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Abstract 

Recent observations of a large excess of cosmic-ray positrons at high energies have 
raised a lot of interest in leptonic decay modes of dark matter particles. Neverthe- 
less, dark matter particles in the Milky Way halo could also decay hadronically, 
producing not only a flux of antiprotons but also a flux of antideuterons. We 
show that for certain choices of parameters the antideuteron flux from dark mat- 
ter decay can be much larger than the purely secondary flux from spallation of 
cosmic rays on the interstellar medium, while the total antiproton flux remains 
consistent with present observations. We show that if the dark matter particle 
is sufficiently light, the antideuteron flux from dark matter decay could even be 
within the reach of planned experiments such as AMS-02 or GAPS. Furthermore, 
we discuss the prospects to observe the antideuteron flux in the near future if 
the steep rise in the positron fraction reported by the PAMELA collaboration is 
interpreted in terms of the decay of dark matter particles. 
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1 Introduction 



There is mounting evidence for the existence of a hitherto undiscovered elementary 
particle in the Universe, the dark matter particle, with an abundance which is approxi- 
mately six times bigger than the abundance of baryons. Dark matter particles are also 
known to exist in galactic halos through their gravitational influence on the motion 
of stars [1]. However, although dark matter particles are ubiquitous in the Universe 
at large, as well as in our own Galaxy, very little is known about their nature and 
properties. 

One of the key questions about the dark matter properties is whether the dark 
matter particles are stable or not. Indeed, astrophysical and cosmological observations 
do not require the dark matter particles to be absolutely stable, but only to have 
a lifetime larger than the age of the Universe. Furthermore, recent evidence for the 
existence of a primary component in the high energy electron and positron fluxes [2] 
could be naturally accounted for by the decay of dark matter particles with a mass 
TTiBM ^ 300 GeV, which decay preferentially into charged leptons of the first and second 
generations with a lifetime Tdm ~ 10^*^s [3]. 

The properties of the dark matter particles are further constrained from observa- 
tions of the diffuse gamma-ray flux by EGRET [4], as well as from observations of 
the antiproton flux by PAMELA [5], BESS95 [6], BESS95/97 [7], CAPRICE94 [8], 
CAPRICE98 [9] and IMAX [10]. More concretely, the good agreement of the theoreti- 
cal predictions for a purely secondary antiproton flux with the measurements indicates 
that the contribution to the total antiproton flux from dark matter matter can only 
be subdominant. On the other hand, the gamma-ray excess in the extragalactic back- 
ground, revealed by the analysis in [11] after subtracting the galactic foreground, hints 
at a possible exotic component of diffuse gamma rays, which will be confirmed or 
refuted by Fermi in the near future. 

In this paper we will estimate the antideuteron flux from dark matter decay as a 
complementary way to determine the properties of the dark matter particle. As pointed 
out in [12,13] , for energies smaller than ~ 3 GeV the expectations for a purely secondary 
antideuteron flux, due to spallation of cosmic rays on the interstellar medium, lie below 
the present BESS hmit on the antideuteron flux [14], as well as the projected limits 
for AMS-02 [15,16] and GAPS [17,18]. Therefore, no detection of antideuterons is to 
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Table 1: Parameters characterizing some commonly used halo models. 



be expected in these experiments if the antideuteron flux is of purely secondary origin. 
At the same time, a discovery of cosmic antideuterons in the upcoming generation of 
experiments would be strong evidence for a primary component and, if interpreted as 
the result of dark matter decay, as evidence for hadronic decays of dark matter particles, 
which may not be inferred from the antiproton flux due to the large astrophysical 
background in this channel. 

In Section 2 we will review antideuteron production from dark matter decay using 
the coalescence model. In Section 3 we will review the propagation of antideuterons 
in the diffusive halo of the Milky Way. In Section 4 we will show our results, and in 
Section 5 we will flnally present our conclusions. 

2 Antideuteron Production 

We will assume that the Milky Way dark matter halo is populated by dark matter 
particles with mass moM, with their distribution following a density profile p(r), where r 
denotes position with respect to the center of the Galaxy. The dark matter distribution 
is usually parametrized as a spherically symmetric profile 



where r = |r| and the parameters a, /3, 7 and Tc are listed in Table 1 for some commonly 
used halo models. Finally, po is a parameter that is adjusted to yield a local halo density 
of p(rQ) = 0.30GeV/cm^ [19], with = 8.5 kpc being the distance of the Sun to the 
Galactic center. 

Dark matter particles at r eventually decay with lifetime Tdm, producing antideuterons 
at a rate per unit energy and unit volume given by 
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where dN^/ dEj; is the energy spectrum of antideuterons produced in the decay. 

The production of antideuterons in the fragmentation of weak gauge bosons or Higgs 
bosons can be described by the coalescence model which assumes that the probability 
in momentum space of producing an antideuteron is proportional to the product of the 
probabilities of producing a single antiproton and a single antineutron: 
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where 7 is a Lorentz factor and C{kp,kn) is the coalescence function which can only 
depend on the relative momentum between antiproton and antineutron, A — kp — kn. 
Using that k^— kp + kn, one obtains 
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The coalescence function is strongly peaked at | A | ~ since the binding energy of 
the antideuteron, B ~ 2.2 MeV, is much smaller than its rest mass, ~ 1.88 GeV. 
Therefore, the previous expression can be approximated by 
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It is useful to define the coalescence momentum, pq, as 
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which can be interpreted as the maximum relative momentum between the antiproton 
and the antineutron for which an antideuteron will form by fusion of the two nucleons. 
The coalescence momentum can be determined experimentally from proton-nucleus 
collisions, yielding po ~ 79 MeV^ [22], or from e+e" collisions at the Z° resonance, 
yielding po ~ 71.8 ± 3.6 MeV [23]. Since the decay of weak gauge bosons is precisely 
the source of antideuterons in the decaying dark matter scenario, we will adopt the 
latter value po ~ 71. 8 ±3.6 MeV in the remainder of the paper. Note that the measured 
coalescence momentum is not far from the estimate obtained by using the antideuteron 



^The uncertainty in the measurement of the hadronic production cross sections translates into a 
range for the coalescence momentum that was estimated to be po = '^^-13 MeV [13]. 
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binding energy, namely y' nipB ~ 46 MeV, supporting the above physical interpretation 
of the coalescence momentum. 

The dark matter decay produces antideuterons, antiprotons and antineutrons with 
an isotropic distribution, and consequently 
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Therefore, Eq. (5) can be cast as: 
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where we have assumed, employing isospin invariance, that the probability of producing 
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an antiproton with momentum kp in the fragmentation is the same as the probability 
of producing an antineutron. 



3 Propagation 

Antideuteron propagation in the Milky Way can be described by a stationary two-zone 
diffusion model with cylindrical boundary conditions [24]. Under this approximation, 
the number density of antideuterons per unit energy, fg{E,f,t), satisfies the following 
transport equation: 

= ^ = V ■[KiE,r)Vh-V,if^fi+QiE,r) . (9) 

where we have neglected non- annihilating interactions of antideuterons with the in- 
terstellar gas, as well as the changes in the energy of the antideuteron during their 
propagation due to energy losses and reacceleration.^ The boundary conditions re- 
quire the solution fd{E, f, t) to vanish at the boundary of the diffusion zone, which is 
approximated by a cylinder with half-height L = 1 — 15 kpc and radius R = 20 kpc. 

The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The 
diffusion coefficient K{E,r) is assumed to be constant throughout the diffusion zone 
and is parametrized by: 

K{E) = /3 7^^ (10) 

^This assumption is well justified for antideuterons from dark matter decay, since after being 
produced in the Milky Way halo, they rarely cross the disk before reaching the Earth. 
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where (5 — v/c and TZ is the antideuteron rigidity, which is defined as the momentum 
in GeV per unit charge, TZ = p{GeV)/Z. The normahzation Kq and the spectral index 
S of the diffusion coefficient are related to the properties of the interstellar medium 
and can be determined from the flux measurements of other cosmic ray species, mainly 
from the Boron-to-Carbon (B/C) ratio [25]. 

The second term is the convection term, which accounts for the drift of charged par- 
ticles away from the disk induced by the Milky Way's Galactic wind. It has axial direc- 

— * 

tion and is also assumed to be constant inside the diffusion region: V^(r) = Vc sign(2;) e^. 
The third term accounts for antimatter annihilation when it interacts with ordinary 
matter in the Galactic disk, which is assumed to be an "infinitely" thin disk with 
half- width h — 100 pc. The annihilation rate, Fann, is given by: 



In this expression it has been assumed that the annihilation cross section between 
an antideuteron and a helium nucleus is related to the annihilation cross section be- 
tween an antideuteron and a proton by the simple geometrical factor 4^/^. On the 
other hand, ~ 1 cm~^ is the number density of Hydrogen nuclei in the Milky Way 
disk, riHe ~ 0.07 the number density of Helium nuclei and cr*j^ is the annihilation 
cross section. No experimental data are available for the dp collisions, although mea- 
surements exist for the total cross section of the charge-conjugated process dp ^ X. 
Since strong interactions preserve charge conjugation, it is reasonable to assume that 
(t|°* = a^J^. However, it is not the total cross section that is required to compute the 
depletion of antideuterons during their propagation, but the annihilation cross sec- 
tion. Unfortunately, there is no experimental information about the annihilation cross 
section neither for dp nor for dp collisions. However, it was noted in [13] that the 
total cross section for dp ^ X can be well approximated, within ~ 10%, by 2(7*^* and 
hence (7|°* ~ '^<^^pp- Then, assuming that the same rule applies for the annihilation 
cross section, one obtains ~ '2.(7^. For our numerical analysis, we will adopt the 
parametrization by Tan and Ng of the proton-antiproton annihilation cross section [26]: 
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36 Tp-o-s mbarn , Tp > 15.5 GeV , 
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where Tp is the kinetic energy of the antiproton. 
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Table 2: Astrophysical parameters compatible with the B/C ratio that yield the minimal 
(MIN), median (MED) and maximal (MAX) flux of antideuterons. 



The solution of the transport equation at the Solar System, r — Tq, z — 0, can be 
formally expressed by the convolution 

1 /-"^DM dNj( E'l) 

ME,) = -—^ / dE'/^,(E,, E'f-^ . (13) 

The analytic expression for the Green's function reads [27]: 

GAT, T) = g exp (-^) ^j^^^-lg^,. (c.|) HT - T') , (14) 



where 



(15) 



and 



A,{T) = 2hT,^^{T) + V, + kS^{T)coth^^^ , (16) 



(17) 



We find that the Green's function can be numerically approximated by the following 
interpolation function: 

Ga{T, T') ~ 10^^ ^x+ylnT+zln-^T^^j,, _ j,^ ^^-3 ^ ^ (^g) 

which is valid for any decaying dark matter particle. The coefficients x, y and z for the 
NFW profile can be found in Table 3 for the various diffusion models in Table 2; the 
dependence of the Green's function on the halo model is fairly weak. In this case the 
approximation is better than about 5% for Tj between 1 and 100 GeV/n and better 
than about 20% for Tj between 0.1 and 1 GeV/n. 
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model 


X y z 


MIN 
MED 
MAX 


-0.3889 0.7532 -0.1788 
1.6023 0.4382 -0.1270 
3.1992 -0.1098 -0.0374 



Table 3: Coefficients of the interpolating function Eq. (18) for the antideuteron Green's 
function for the NFW halo profile. 

Finally, the flux of primary antideuterons at the Solar System from dark matter 
decay is given by: 

^f^iEs) = l^h{E,l (19) 

where v is the velocity of the antideuteron. 

However, this is not the antideuteron flux measured by balloon or satellite exper- 
iments, which is affected by solar modulation. In the force field approximation [28] 
the effect of solar modulation can be included by applying the following simple for- 
mula that relates the antideuteron fiux at the top of the Earth's atmosphere and the 
interstellar antideuteron fiux [29]: 

^TOA/^ \ _ ( 2m jTtoa + T^oA \ ^is/^ \ /oan 
2m,T:s + ri )^^^^'^^^ (20) 

where Tis = Ttoa + ^f, with Tjs and Ttoa being the antideuteron kinetic energies at 
the heliospheric boundary and at the top of the Earth's atmosphere, respectively, and 
0ir being the solar modulation parameter, which varies between 500 MV and 1.3 GV 
over the eleven-year solar cycle. Since experiments are usually undertaken near solar 
minimum activity, we will choose 0f = 500 MV for our numerical analysis in order to 
compare our predicted flux with the collected data. 

4 Antideuteron Flux at Earth from Dark Matter 
Decay 

Using the injection spectrum of antideuterons from dark matter decay calculated in 

Section 2 and the propagation formalism described in Section 3, it is straightforward to 

calculate the antideuteron flux at Earth from dark matter decay. ^ We will pursue here 

■^For calculations of the antideuteron flux from the annihilation of dark matter particles, see 
[12,13,30]. 
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a model-independent approach, calculating the expected antideuteron flux at Earth for 
various hadronic decay channels and different dark matter masses. To better compare 
the predictions for the different possibilities, we will first fix the dark matter lifetime 
to be tdm = lO^^s, which is the order of magnitude which could explain the HEAT 
and PAMELA anomalies in the positron fraction [31,32,3,33] and which saturates the 
EGRET constraints on the diffuse gamma-ray flux [34,32]. Later on, we will discuss 
in more detail the predictions for the antideuteron flux in the light of the PAMELA 
anomaly. 

In the case that the dark matter particle is a fermion ■0, the following hadronic 
decay channels are possible: 

^ . (21) 

The fragmentation of the weak gauge bosons and the Standard Model Higgs boson 
produces a flux of antideuterons which has been calculated using the coalescence model 
to simulate the nuclear fusion of an antiproton and an antineutron (c/. Section 2), 
and the event generator PYTHIA 6.4 [35], to simulate the production of individual 
antiprotons and antineutrons. 

Clearly, an exotic contribution to the antideuteron flux is correlated with an ex- 
otic contribution to the antiproton flux, which is severely constrained by a number 
of experiments. Namely, the measurements of the antiproton flux by PAMELA [5], 
BESS95 [6], BESS95/97 [7], CAPR1CE94 [8], CAPR1CE98 [9] and IMAX [10] do not 
show any signiflcant deviation from the predictions by conventional astrophysical mod- 
els of spallation of cosmic rays on the Milky Way disk. Therefore, in order to evaluate 
the prospects of observing the antideuteron flux in future experiments, one has to en- 
sure that the prediction for the total antiproton flux does not exceed the observed flux. 
The antiproton flux from dark matter decay has been calculated following the analysis 
in [31], for the same propagation models as for the antideuteron flux. 

We show in Figs. 1, 2, 3 the antiproton and antideuteron fluxes from dark matter 
decay as a function of the kinetic energy per nuclcon, assuming that the dark matter 
particle is a fermion ip, which decays exclusively into Z^u, W^i^ or h^u, respectively. 
The lifetime has been flxed to tdm = lO^^s, while the dark matter mass has been 



8 



chosen to be uibm — 200, 400, 600, 800 GeV. The predictions for the antiproton and 
antideuteron fluxes for a general model with a different dark matter lifetime and with 
arbitrary branching ratios can be straightforwardly derived from these figures. In each 
plot, wc show the prediction for the antiproton and antideuteron fiuxes at the top of 
the atmosphere for the MIN, MED and MAX propagation models (see table 2) as well 
as the expected flux from spallation, taken from [27] in the case of the antiprotons 
and from [13] in the case of the antideuterons.^ For the halo model, we adopted the 
Navarro-Prenk- White (NFW) proflle; the results for other halo proflles are very similar 
to the ones presented here. We also show the present upper limit on the antideuteron 
flux from BESS [14], as well as the projected hmits for AMS-02 for three years of data 
taking [15,16], as well as for GAPS [17,18] for a long duration balloon (LDB) flight 
(60 days total over three flights) and for an ultra-long duration balloon (ULDB) flight 
(300 days total). 

For this particular choice of the lifetime, t^m — 10^^ s, we find wide ranges of 
diffusion parameters yielding a total antiproton fiux consistent with the observations. 
Besides, the total antideuteron flux always hes well below the present BESS bound. We 
also flnd that for dark matter masses below ~ 1 TeV, the antideuteron flux from dark 
matter decay can be, at energies below 3 GeV, signiflcantly larger than the secondary 
antideuteron flux from spallation. Furthermore, the primary antideuteron flux from 
dark matter decay could be large enough to be observable at the projected experiments 
AMS-02 and GAPS. Therefore, since the purely secondary antideuteron flux is expected 
to be below the sensitivity of projected experiments, the observation of an antideuteron 
flux in the near future could be interpreted as a signature of dark matter particles which 
decay hadronically. 

When the decaying dark matter particle is a scalar, the following hadronic decay 
modes are possible: 

^ h°h°. (22) 

*In these plots we only show the prediction for the secondary antiproton and antideuteron fluxes 

for the MED propagation model. Whereas the uncertainty on the propagation model is typically not 
very important for the prediction of the secondary fluxes, there exists a more important source of 
uncertainty stemming from the nuclear and hadronic cross sections, which can be as large as 25% for 
antiprotons [27] and 100% for antideuterons [13]. 
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Figure 1: Fluxes of antiprotons (upper curves) and antideuterons (lower curves) from the 
decay of fermionic dark matter particles, assuming that the dark matter particle decays 
exclusively as ip ^ Z^v with a fixed lifetime tdm = lO'^^s. The dark matter mass is 200 
GeV (top left panel), 400 GeV (top right), 600 GeV (bottom left) and 800 GeV (bottom 
right). The dashed lines indicate the primary fluxes from dark matter decay for the MIN 
propagation model, the dotted lines for the MED propagation model, and the dash-dotted 
lines for the MAX propagation model (c/. Table 2). The solid lines indicate the secondary 
fluxes. 

In these cases, the produced antideuteron flux is approximately twice as large as in 
the fermionic channels, — Z^v^ W^t^ ^ ifj h^u, respectively, and will not be 
discussed further. 

Let us now discuss the prospects for the detection of an antideuteron flux in pro- 
jected experiments, assuming that dark matter decay is the explanation of the positron 
excess reported by the PAMELA and HEAT collaborations. As discussed in [3], the 
steep rise of the positron fraction measured by PAMELA can be explained by the decay 
of a dark matter particle into hard electrons and positrons. Namely, in the case that the 
dark matter particle is a fermion, the decay modes ip —>■ e^e~v and ip PF^e^ (and 
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Figure 2: The same as Fig. 1 for the decay mode ifj W^i"^. 



the analogous decay modes into muon flavor) are favored by the data. On the other 
hand, when the dark matter particle is a scalar, the decay modes (p e~^e~ , fi'^fj^' are 
favored. Purely leptonic decay modes, such as ip ^ or — > do not produce 

antideuterons. However, as discussed above, the decay of a fermionic dark matter par- 
ticle into bosons and charged leptons could produce an observable antideuteron 
flux. 

We show in Fig. 4 the total antiproton and antideuteron fluxes from the decay of a 
fermionic dark matter particle in the channel ip — > W^e^, flxing for each dark matter 
mass the lifetime in order to account for the steep rise in the positron fraction observed 
by PAMELA. For a dark matter mass ttidm = 300, 600, 1000 GeV, the corresponding 
lifetimes are tdm = 4.0, 2.3, 1.6 x 10^^ s, respectively. For these particularly interesting 
decaying dark matter scenarios, we flnd that the antideuteron flux could be within 
the reach of the planned experiments AMS-02 and GAPS, provided the dark matter 
particle is not too heavy. The discovery of antideuterons would thus favor the decay 
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Figure 3: The same as Fig. 1 for the decay mode t(j hPv. 

mode if) — > W^e^ as a possible origin of the PAMELA positron excess over the purely 
leptonic decay modes such as ip ^ e^e~v.^ It is important to emphasize that this 
conclusion holds for ranges of propagation parameters which not only reproduce the 
positron excess observed by PAMELA, but also yield a total antiproton flux consistent 
with present measurements. 

5 Conclusions 

We have calculated the antideuteron fluxes at Earth from the decay of dark matter 
particles in the Milky Way halo and we have discussed the prospects to observe cos- 
mic antideuterons in the projected experiments AMS-02 and GAPS. We have adopted 
a model-independent approach analyzing possible hadronic decay modes for the dark 
matter particle. The nuclear fusion of an antiproton and an antineutron to form an- 
^The decay mode tp ~^ W^e^ may also yield signatures in the diffuse gamma-ray background. 
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Figure 4: Total antiproton and antideuteron fluxes including a primary contribution to 
the flux from the decay of a fermionic dark matter particle in the channel '0 W^e^ . The 
lifetime has been chosen, for each dark matter mass, to reproduce the steep rise in the positron 
fraction observed by the PAMELA collaboration. When the dark matter mass is 300 GeV 
(top left panel), 600 GeV (top right panel) and 1000 GeV (bottom panel), the corresponding 
lifetimes are tdm = 4.0, 2.3, 1.6 x 10^^, respectively. The dashed lines indicate the total 
fluxes for the MIN propagation model, the dotted lines for the MED propagation model, and 
the dash-dotted lines for the MAX propagation model (c/. Table 2). We also show the purely 
secondary fluxes as solid lines. 
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tideuteron was simulated employing the coalescence model, while the propagation of 
antideuterons in the Galaxy was described by a stationary two-zone diffusion model 
with cyhndrical boundary conditions. The hadronic showers also produce a flux of pri- 
mary antiprotons, which is severely constrained by observations. Therefore, wc have 
simultaneously calculated the predicted antiproton flux, to verify whether the chosen 
parameters and decay modes are consistent with present observations. 

We have shown that there are choices of parameters where the antideuteron flux 
from dark matter decay can be much larger than the purely secondary component 
from spallation of cosmic rays on the interstellar medium, while at the same time 
the total antiproton flux remains consistent with observations. We have also shown 
that if the dark matter particle is sufficiently light, the antideuteron flux from dark 
matter decay could be within the reach of the planned experiments AMS-02 or GAPS, 
while the secondary component is expected to lie below the sensitivity of planned 
experiments. Therefore, the observation of cosmic antideuterons in the near future 
could be interpreted as an indication for hadronic decays of dark matter particles. 
In particular, this conclusion holds for a fermionic dark matter particle which decays 
preferentially via — > W^e^, which has been proposed as an explanation for the steep 
rise in the positron fraction observed by the PAMELA collaboration. 
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